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The bag parameters and the decay constants of neutral Bt s \ and K mesons were among the non- 
perturbative hadronic inputs to the classical CKM Unitarity Triangle Analysis. Thanks to the big 
amount of experimental information collected in the last few years at the B-factories and by the 
CDF collaboration, these matrix elements are now among the outputs of the unitarity fits, once the 
validity of the Standard Model has been postulated. Lattice calculations of the mixing amplitudes 
are still needed in order to make a test of the theory, provided that their statistical and systematic 
errors are under control at the level of a few percent. Here we review some of the recent lattice 
calculations of these quantities. 



I. INTRODUCTION 



Lattice QCD calculations of the B/ s \ and K mixing 
amplitudes where needed in the past in order to be used 
as inputs to the Unitarity Triangle Analysis (UTA). Af- 
ter the recent measurements of the angles of the unitarity 
triangle the p-r\ plane is over-constrained and the mixing 
amplitudes can now be extracted by one of the possi- 
ble variants of the UTA (see [l[ and references there). 
From the phenomenological point of view it is thus legit- 
imate to ask whether it is still needed a precise lattice 
calculation of the mixing amplitudes. In answering this 
question the point to note is that the unitarity fits are 
over-constrained only within the Standard Model. In this 
scenario one can make a precision test of the Standard 
Model and hopefully reveal the presence of new physics 
by comparing the "experimental" determinations of the 
mixing amplitudes with their theoretical predictions. To 
this end lattice is still needed but it has to provide, let's 
say, AM S with an error of the same order of magnitude 
of the one quoted by the CDF collaboration, i.e. all the 
systematic errors have to be under control at a level of a 
few percent. Among all the systematics that presently af- 
fect lattice calculations 1 the worst are the "uncontrolled" 
ones, in the sense that cannot be reliably quantified with- 
out performing independent calculations with different 
discretizations of the continuum action. To this category 
certainly belong quenching and, in our opinion, rooted 
staggered fermions. Staggered fermions are introduced 
on the lattice by simulating a quark action that suffers 
from doubling, i.e. it has 16 one-component fermions in 
the classical continuum limit that are packaged into 4 
tastes of 4-component Dirac fermions. Rooting means 
that gauge configurations are generated by taking the 
fourth root of the staggered quark determinant 2 . If 



'Electronic address: nazario.tantalo@roma2.infn.it 

1 see C. Sachrajda and T. Onogi talks at this workshop for an 
extended discussion on the lattice systematics 

2 This formalism has been largely used because it is particularly 
cheap from the computational point of view but recently it has 
been proved the feasibility of large scale Nf = 2 lattice simula- 



"taste symmetry" is exact, like in the continuum theory 
rooting is obviously a kosher operation. At finite lattice 
spacing taste symmetry is broken and rooting breaks lo- 
cality and unitarity. Whether these properties are recov- 
ered or not by taking the continuum limit is the subject 
of a strong debate within the lattice community. At the 
annual lattice conference S. Sharpe has given a beautiful 
critical review on the subject Q entitled "Rooted stag- 
gered fermions: Good, bad or ugly?". His answer to this 
question is "ugly, in the sense that they are affected by 
unphysical contributions at regulated stage that need a 
complicate analysis to be removed" . We share his view. 



II. B A 

The so-called bag parameter Bk parametrizes the mix- 
ing amplitude of the Kq-Kq mesons according to 

(K^K ) = (K°\d vv+AA \K°) = *-M 2 K f K B K {jfi 

where the operator 0\ = s*7 M (l — 7s)rf* ^7^(1 — 75)^ is 
usually conveniently decomposed into a parity-even and 
a parity-odd part, 0\ = Ovv+aa — Ova+av- 

The breaking of chiral symmetry on the lattice with 
Wilson fermions complicates considerably the renormal- 
ization pattern of the Ovv+aa operator that happens 
to mix with the other 4 operators in the parity-even ba- 
sis [f| 0, @]. The issue in such a calculation consists 
in obtaining the numerical accuracy required in order to 
keep under control the mixing subtractions or in devising 
smart strategies to avoid the problem. Two groups have 
been able to circumvent the mixing problem. The au- 
thors of refs. [§, have used the so called twisted mass 
formulation of lattice QCD (tmQCD) in order to map the 
matrix element of the parity even operator Ovv+aa into 
the matrix element of the parity odd operator Ova+av 
that renormalizes multiplicatively; they have been able to 



tions with pions as light as 300 MeV and physical volumes of the 
order of 2 fm. 
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FIG. 1: Top. Average of the uncorrelated Bk results without 
the quoted errors. Bottom. Average of the Bk results with 
Nf > with the quoted errors (red lines). 



calculate Bk by keeping under control all the sources of 
errors apart from quenching (non-perturbative renormal- 
ization, estimate of ST/(3)-breaking effects, continuum 
limit with 5 lattice spacings, estimate of finite volume ef- 
fects). The authors of refs. [l(| [H| avoided the mixing by 
using a chiral ward identity that again relates the matrix 
elements of Ovv+aa to that of Ova+av at the price of 
computing on the lattice a four-point Green function. 

In the case of lattice discretizations that satisfy the 
so called Ginsparg- Wilson (GW) relation an exact chi- 
ral symmetry is preserved also at finite lattice spacing. 
Domain wall fermions satisfy the GW in the limit of an 
infinite fifth dimension. Practically, the fifth dimension 
is finite and the lattice chiral symmetry is only approxi- 
mately preserved. The authors of refs. [Ill, have per- 
formed a calculation of Bk with respectively Nf = 2 and 
Nf = 2+1 flavours of dynamical domain wall fermions. 
The Nf = 2 + 1 results have been obtained at fixed lattice 
spacing (a ~ 0.12 fm), with non perturbative renormal- 
ization (by neglecting the small mixing due to the "resid- 
ual mass term"), by interpolating the physical K meson 
state, on a single volume (L ~ 2 fm) ; a simulation at the 
same lattice spacing on a larger volume (L ~ 3 fm) is 
under way. 

On the one hand, there have been so many different 
calculations of Bk among the years that it is not possi- 
ble to enter into the details of all of them in this short 
review 3 (see TABLE QJ. On the other hand none of this 
calculations is able to take under control all the sources 



we have just mentioned some representative calculations and 
apologize with the authors whose results have not been covered 
in greater detail. The same holds also for the following sections. 



TABLE I: Lattice calculations of the renormalization group 
invariant (RGI) kaon bag parameter Bk- 



collaboration 


B K 


N f 


JLQCD97 [12] 


0.868(59) 





BecirevicOO [20] 


1.01(9) 





CP-PACS01 [13] 


0.795(29) 





SPQCDR02 [10] 


0.91(9) 





BosMar03 [14] 


0.87(8) 





MILC03 [15] 


0.79(9) 





BabichOf! llfil 


0.79(8) 


o 


ALPHA06 [18] 


0.735(71) 





RBC03 [21] 


0.697(33) 


2 


UKQCD04 [19] 


0.67(18) 


2 


SPQCDR05 [11] 


1.02(25) 


2 


RBC05 \17} 


0.78(7) 


2 


RBC-UKQCD06 [22] 


0.778(36) 


2+1 


HPQCD-UKQCD06 [23] 


0.85(12) 


2+1 



of systematics at the same time. Since different numbers 
have been obtained with different actions, techniques, as- 
sumptions, etc. we can get an estimate of the systemat- 
ics by averaging all the results that are "uncorrelated" 
(in the sense that we neglect results that have been up- 
dated by the same collaboration at fixed Nf) without the 
quoted errors (see FIG. [T]top plot). As a result we get 
Bk = 0.81(3) i.e. a relative error of the order of 4%; if 
instead we take the average of the numbers with Nt > 
by trusting the quoted errors we get Bk = 0.78(2) (see 
FIG. [T] bottom plot). The previous numbers have to be 
taken as "provocative" averages: unless a clear statement 
is made on which lattice results can be trusted and which 
have to be excluded from phenomenological analysis one 
should conclude that Bk is presently predicted by the 
lattice with a few percent error. A conservative estimate 
of the errors, to be used in phenomenological applica- 
tions, can be obtained for example by accounting for the 
dispersion of the results: 

B K = 0.78(2)(9) (1) 
HI. f Bq 

The decay constants of the B q mesons, where q stays 
for either a down or a strange quark, enter in the 
parametrization of the B q -B q mixing amplitudes together 
with the bag parameters, 

(B q \6 vv+ AA\B q ) = |m|,/| s £b» 

What it is actually needed in order to perform the UTA 
is the combination /r \/B~b 9 > that comes out to have a 
smaller statistical error on the lattice w.r.t. the product 
of fg and \JBs q computed separately. Since there are 
many more calculations of the decay constants than the 
bag parameters and since we want to use as much infor- 
mation as possible in taking the averages, we will discuss 
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separately the lattice calculation of fs , in this section, 
and of Bg q , in the next section, but first we briefly com- 
ment on the issues related to the simulation of the heavy 
flavours on the lattice. 

On currently affordable lattice sizes (at least in un- 
quenched simulations) one has ami > 1 and Lm q > I 
or amj < I and Lm q < I , i.e. a relativistic beauty- light 
meson can be simulated on big volumes with big cutoff 
effect or on small volumes with big finite volume effects. 
The different approaches that have been devised in order 
to solve this problem can be divided into "big volume" 
and "small volume" strategies. 

Big volume strategies. Lattice simulations of heavy 
quarks on physical volumes are performed by recurring 
to effective field theories. One possibility is to simulate 
the lattice static action and, eventually, the relativistic 
theory with heavy quark masses in the charm region in 
order to interpolate the bottom region. Another pos- 
sibility is the so called Fermilab approach [24], [25|, [26[ 
that consists in improving the heavy quark action with 
mass dependent coefficients, provided that \ap\ <C f ; the 
procedure smoothly interpolates between heavy and light 
quarks and the continuum limit can be taken although 
the mass dependence of the improving coefficients makes 
the procedure highly non trivial. Still another possibility 
is to simulate on the lattice a discretized form of the non 
relativistic heavy quark action expanded to a given order 
in v 2 and a s (see for example [271]); as a consequence of 
the non-renormalizability of the theory, the continuum 
limit cannot be taken and the matching with full QCD 
can be done only perturbatively; furthermore the expan- 
sion is expected to work for onium systems but it has 
been widely used also to study heavy-light mesons. 

Small volume strategies. The step scaling method 
(SSM) has been proposed in [28[ in order to deal with 
two-scale, Ei <C Eh, problems in lattice QCD. The 
method starts from the following identity 



0{E h ,E h oo) = 

O(E h ,E u 2L ) O(E h ,E h 4L ) 
O(E h , El ,L )x 0{EhjEhLo) * 0{EhiEli2Lo) - 

* v ' V v ' 

<r(E h ,E,,L ) a(E h ,Et,2L ) 



The idea is to start the computation on a small volume 
where the high energy scale can be set to its physical 
value and to correct for the finite volume effects by com- 
puting the step scaling functions, a(Eh, Ei, L), stopping 
the recursion when the last factor is equal to one within 
the required precision. The strength of the method can 
be better understood by specializing the previous iden- 
tity to the case of a heavy-light meson observable, say 
fg . In this case the dependence of the step scaling func- 
tions upon the b quark mass can be predicted by using 



Fermllab 97 
MILC 98 
JLQCD 99 
UKQCD 00 
APE 00 
ALPHA 03 
ROMEII 03 
RQMEII-ALI 

CP-PACS 
CP-PACS 
MILC 02 
JLQCD 03 
UKQCD 04 
Gsriiysk iji-, 
ALPHA 06 

HPQCD 05 




• 






£ * 


« 8 - 


* e e 






Fermllab 97 q 
MILC 98 
JLQCD 99 
UKQCD 00 
APE 00 

ALPHA 03 ui 

ROMEII 03 

ROM Ell -ALPHA 06 

CP-PACS 00 
CP-PACS 01 5 
MILC 02 
JLQCD 03 
UKQCD 04 
Gadiyak 05 

ALPHA 06 J 
HPQCD 05 




«. 

«. 






t) c « 




: r— ? 

" i i * * * » « 


Mi 





FIG. 2: Top. Average of the uncorrelated /s 3 results without 
the quoted errors. Bottom. Average of the /s s results with 
Nf > with the quoted errors (red lines). 



HQET both at numerator and denominator 
f° B (m d ,2L) (l 



a(m h ,m d ,L) 



/i (m d ,2L) 



stat 



(m d ,L) 1 + 



f B (m d ,2L)-f B {m d ,L) , 



m h 



The maximum number of points of the lattice discretiza- 
tion is fixed during the recursion so that, in order to 
have arrih < 1 at each step, the step scaling functions 
have to be computed at heavy quark masses smaller than 
the physical beauty mass. Nevertheless, from the previ- 
ous equation, it is clear that the step scaling functions 
depend mildly upon the high energy scale thank to the 
cancellation f B (m d ,2L) — f B (m d ,L) that gets stronger 
and stronger as the volume is increased. Another use- 
ful application of finite volume techniques has been used, 
originally in ref. (29| to implement the renormalization 
procedure of HQET fully non-perturbatively. The idea 
is to match the effective theory with full QCD on a small 
volume at the scale /x = m&. Then the running of the 
matching factors is computed through a step scaling re- 
cursion in the effective theory. This method can be ap- 
plied in conjunction with the SSM to compute, for exam- 
ple, cr stat (m d ,L) (see ref. 0). 

We now come to the lattice results for f Bs and 
r_e = frt -//«■ Among the quenched results those of 
refs. [30, [42] have been obtained through the SSM by 
keeping under control all the systematics apart from the 
quenching. Among the Nf — 2 results, the ones of 
refs. [iH l5u| have been obtained within the static approx- 
imation^ ving a particularly high value of r B ; also the re- 
sult of [52| is static but non-perturbatively renormalized. 
The Nf = 3 results [i^, [53[ have been obtained by using 
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TABLE II: Lattice calculations of f Bs (MeV) and r B = 
f B Jf B . 



collaboration 


}b s 




Nf 


Fermilab97 [38] 


185(16) 







MILC98 [39] 


171(44) 







JLQCD99 [40] 


191(18) 







UKQCD00 [32] 


204(27) 







APE00 [33] 


235(20) 







ALPHA03 [45] 


206(10) 







ROMEII03 [47] 


192(7) 







ROMEII-ALPHA06 [30] 


191(6) 







CP-PACS00 [41] 


250(18) 


1.203(64) 


2 


CP-PACS01 [42] 


242(52) 


1.179(29) 


2 


MILC02 [43] 


217(36) 


1.16(5) 


2 


JLQCD03 [36] 


215(17) 


1.13(12) 


2 


UKQCD04 [48] 


256(45) 


1.38(15) 


2 


Gadiyak05 [50] 


341(32) 


1.38(15) 


2 


ALPHA06 [52] 


297(14) 




2 


HPQCD05 [49] 


259(32) 




2+1 


Fermilab-MILC-HPQCD06 [53] 


253(42) 


1.27(6) 


2+1 



the mixing, as for Bk, and the problem related to the 
presence of a heavy and a light quark, as for f Bq . For 
this reason the number of lattice calculations of B Bq is 
much smaller than in the case of f Bq or Bk- 

Nevertheless, by looking at TABLE lTTTl it emerges that, 
within the quoted errors, Bb does not seem to depend 
upon the number of dynamical flavours, the rcnormal- 
ization systematics (the quenched result of ref. [34[ has 
been non-perturbatively renormalized), the strategy used 
to handle with heavy quarks and even the light quark 
mass. Actually, the matrix element {B q \Ovv+AA\B q ) 
does show a sizable dependence upon all these variables 
but through the vacuum saturation approximation, i.e. 
8Af| 5 /^/3. The average of the Nf > calculations, 
with an error that takes into account the spread of the 
results, is 

B Ba (m b ) = 0.84(3)(5) B B (m b ) = 0.83(1)(6) (3) 



V. A CALCULATION OF G(w) 



TABLE III: Lattice calculations of BB s (m,b) and BB{rn,b) 



collaboration 


BB s (m b ) 


B B (m b ) 


N f 


UKQCD00 [32] 


0.90(4) 


0.91(6) 





APE00 [33] 


0.92(7) 


0.93(10) 





SPQCDR01 [34] 


0.87(5) 


0.87(6) 





JLQCD02 


35] 


0.86(5) 


0.84(6) 





JLQCD03 


36] 


0.850(64) 


0.836(68) 


2 


Gadiyak05 [50] 


0.864(76) 


0.812(82) 


2 


HPQCD06 [37] 


0.76(11) 




3 



rooted staggered fermions for the dynamical light quarks 
and NRQCD and Fermilab respectively for the heavy. By 
looking at FIG [2] it emerges that quenched lattice calcu- 
lations, though compatible within themselves, are sys- 
tematically smaller than unquenched results; the quoted 
errors are still large but unqucnching seems to have a sig- 
nificant effect on this observable (this is not the case of 
Bk within the quoted errors) and the "provocative" av- 
erage (f Bs = 245(13) MeV and r B = 1.24(4)) prefers un- 
quenched results. The unquenched average, static points 
included, with an error that takes into account the spread 
of the results is 



f Ba = 268(17)(20) MeV, 



fB 




1.20(2)(5) (2) 



FIG. 3: Comparison of lattice data and experimental deter- 
minations of G(lo): Vet = 41(5) x 10~ 3 is extracted by using 
the experimental points at uj = 1.2. 

We now change subject to put up to the results of a 
preliminary quenched calculation [54] of the form factor 
G(uj), needed in order to extract V c b from the exclusive 
scmileptonic decay B^ — > D^iv (uj = p B -p B / M B Mu)- 
The calculation has been carried on by using the SSM and 
by defining the form factor and the kinematical factors in 
terms of ratios of three point correlation functions. The 
results have been obtained with a relative error of about 
4% for values of 1 < u> < 1.2 were experimental data do 
not need to be extrapolated (FIG. [3]). 
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